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ABSTRACT 


Motivated by the recent explosive development of results in the area of parametric 
robust control, this paper presents a new technique to identify a family of uncertain 
systems. The new technique takes the frequency domain input and output data 
obtained from experimental test signals and produces an “interval transfer function” 
that contains the complete frequency domain behavior with respect to the test signals. 
This interval transfer function is one of the key concepts in the parametric robust 
control approach and identification with such an interval model allows one to predict 
the worst case performance and stability margins using recent results on interval 
systems. The algorithm is illustrated by applying it to an 18 bay Mini-Mast truss 
structure. 

1. INTRODUCTION 

It is well known that obtaining an accurate mathematical description of a system is 
impossible, usually very costly and often increases the complexity of the corresponding 
control mechanism. Robust control deals with systems described as a family consisting 
of a nominal model with uncertainty around it. In particular, parametric robust 
control deals with systems whose parameters of interest vary in known independent 
intervals. The simplest form of such a system is called an “interval system” wherein 
the model consists of a family of transfer functions with numerator and denominator 
coefficients varying in prescribed independent intervals. 

A recent trend in the area of system identification is to move toward a new 
direction in which one tries to model the system uncertainties to fit the available 
analysis and design tools of robust control. The motivation for this direction of 
research is clearly stated in [1]. The excellent collection of papers dealing with this 
important issue extensively is also found in [2] and references therein (also see [3,4,5]). 
The main emphasis of these results is either to classify the plant uncertainty as a 
various form of norm bounded uncertainties such as H°° } /i, and structured singular 
values or to find the optimal nominal models to fit their respective robust control 
design methods. 

Commonly, an interval transfer function is interpreted as a family of transfer 
functions whose coefficients are bounded by some known intervals and centered at 
the nominal values. In many cases, this type of framework is unnatural and physical 
parameter perturbations do not correspond to transfer function coefficients. In order 
to relax this limitation, approaches to deal with linearly or multilinearly correlated 
perturbations have been proposed. On the other hand, if we observe recent devel- 
opments in the interval system area we see that the nominal system has very little 
significance. These results in fact emphasize the boundary properties of the family 
of systems under consideration. In fact, virtually all the important results we enjoy 
today in this field are based on the boundary generating extreme points, edges, and 
segments of the interval system. 

Following the motivation described in [1], suppose that the behavior of the plant 
is described by some known test input and its corresponding measurement output. 
Due to many reasons such as noise, inaccurate measurements etc., a fixed linear time- 
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invariant identified model will almost never exactly represent the data obtained from 
the plant. Our aim in this paper is to obtain a reasonable interval transfer function 
model around (not necessarily centered in) the nominally identified transfer function 
so that the entire frequency domain behavior of the physical plant is completely con- 
tained in that of the model. We give a systematic and simple algorithm to accomplish 
this task. We also discuss the issue of design validation of the obtained interval trans- 
fer function. Finally we demonstrate this practical technique by applying it to an 18 
bay Mini-mast truss which has been extensively used for various large space structure 
control problems. 

2. PROBLEM FORMULATION 

Consider the following configuration shown in Figure 1. 


test input 


system 


Figure 1. 


test output 


In system identification one applies test inputs and measures the system response 
in order to identify the parameters of an appropriate proposed mathematical model. 
It is also common that these test signals are represented in a form of frequency domain 
data. In system identification literature there are numerous techniques available (see 
[7] and references therein) to determine a best possible linear time - invariant model 
that fits this set of frequency domain data. Suppose that the test frequencies are 
■ '&N and the complex numbers u(juji), y(jwi) denote in phasor notation the 
input-output pair at the frequency u Let 

y(M) •= i = 1, 2, • • ■ , N (2.1) 

denote the test data generated from an identification experiment. Suppose that G T {s) 
is the transfer function of a linear time-invariant system which is such that G J (ju;) 
is closest to D(ju ;) in some norm sense. In general it is not possible to find a single 
rational function G T (s) for which G z ( j&i) = D(jtOi) and the more realistic identifica- 
tion problem is to in fact identify an entire family G(s) of transfer functions which is 
capable of “explaining” the data in the sense that for each data point there 

exists some transfer function Gi(s) 6 G(«s) with the property that Gi(ju>i) = 

The family G(«s) can be parametrized in many alternative ways. For instance an un- 
structured approach to describing G(,s) using a normed algebra is to let each element 
G(s) of G(s) be described as G(s) = G I (s) + A G(s) where the norm | AG($)| < p . In 
such a case the family G(s) is identified once G^(s) and p are determined. In general 
the identification algorithm should also be efficient in the sense that the family G(s) 
that it produces should be ideally minimal among the set of all such families that 
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explain the data. In the unstructured case described above this translates to choosing 
a small value of p. 

In this paper our objective is to develop an identification algorithm in a frame- 
work where the family of linear time-invariant systems G (s) is obtained by letting 
the transfer function coefficients lie in intervals around (not necessarily centered in) 
those of the nominal G I (s ). Of course the identification requirement is that 

D(jui) G G (juji) V (2.2) 


Let 

G'W 

Uq + UiS + 7i2<S 2 + “1" * * 

■ + n n s n 

(2.3) 


d 0 + d-is + d 2 s 2 + d 3 s 3 + • • 

• + d n s n ' 

We define 

GW 

h 0 + his + h 2 s 2 + h 3 s 3 + ■ • 

■ + h n s n 

(2.4) 


do + d\S + + d^S^ + • ■ 

• + d n s n 

and 

GW ~ {GW : 

where 

: hi € [rii 

- + w ni e+],<fi E [d { - 

W * e *> rf + ;]> 

V t} 
(2.5) 


w := 

W do ■ 

■ - W dn 

^no 

W n„ 

e + := 

et ■ 
“° 

‘ e d 

“n 

6 no 

• ■ ' < ] 

:= 

. e do • ■ 

” e d n 

£ no 

... e - 

Tl n J 


The components of w are to be regarded as weights chosen apriori whereas the es are 
to be regarded as dilation parameters to be determined by the identification algorithm 
and the data D(ju>i) 

Remark 1. Note that in the expression in eq. (2.5) we use vectors and e^. instead 
of a single e. This setting allows and d± to not necessarily be the center point of the 
intervals in which hi and d { , lie respectively. This flexibility is important to achieve 
the minimum possible size of the family G(s). 

The requirements on the identified interval model G(s) become: 

1) Membership Requirement: D(juii) E G(jf’u^) V x. 

2) Size Requirement: ||e ± || as small as possible. 

3) Frequency Response Requirement: the weights must be chosen so that the 

frequency response of G(ju>) is bounded as tight as possible for every frequency. 

It is important to note that both size and frequency response requirements are 
crucial because smaller intervals do not necessarily map to smaller image sets. The 
frequency response requirement ensures that the frequency domain image set be as 
small as possible. 
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3. INTERVAL SYSTEM MODELLING 


Our objective is to find the set, an interval system G(s), to satisfy the three require- 
ments given in the previous section. As described above the procedure is divided 
into two part. First, we identify a linear time - invariant model G T (s) which repre- 
sents the test data D(ju>) as closely as possible. A variety of algorithms are available 
in the system identification literatures. We use a least squares algorithm which is 
widely known [7]. Then using this identified model as a nominal model, we create 
the tightest intervals around each coefficient of the nominal transfer function G J (s) 
while satisfying the membership and frequency response requirements. 


3.1. Nominal System Identification 

In this subsection, Since the purpose of this paper is not to deal with traditional single 
system identification problem, we briefly describe a standard method to identify a 
nominal transfer function whose frequency response fits the given test data 
as closely as possible. The method we use here is widely known as the weighted least 
square method [7]. The least square approach to curve fitting a transfer function in 
the frequency domain may be found in [8]. An appropriate order of model may be 
determined by checking the singular values of the Hankel matrix generated from the 
impulse response data. Under the assumption that the data is noise free, the number 
of nonzero singular values determines the order of the system. The details of this 
approach is found in [9]. After determining the appropriate order of the system, we 
let the nominal transfer function be 




n ( s ) 

d(s) 


(3.1) 


The nominal transfer function coefficients must be selected to minimize the following 
index: 


E - n(M))} 2 + {Im[D(M)d(M) - niM)}} 2 } . (3.2) 

i— 1 

This least square problem generates 2N linear equation for 2 n unknown coefficients 
of the transfer function. The weight W 1 may be selected by finding the minimum 
variance estimator of unknowns. The details of methods for selecting the normalized 
least square weight W 1 {jw) and its effect on the identified model are discussed in [7]. 
The reference [7] also discusses the problem of selecting an appropriate order of the 
transfer function in detail. In general, the relative error in the valley parts of the 
frequency response is more significant than the one in the peak parts, it is necessary 
to assign high weights for the frequency ranges in the valley parts of the frequency 
response. 

Using this idea, we select the normalized least square weight W T (ju) and improve 
the identified model by comparing the identified model and the test data. 
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3.2. Weight Selection 

As shown in eq. (2.5), the size of the interval of variation for each coefficient of the 
family G(s) depends on w and e. In this subsection, we consider the problem of 
finding an appropriate set of weights w. The weight selection procedure is extremely 
important because inappropriate selection of weights may result in an unnecessarily 
large family. This results in a large image set in the complex plane, even though the 
intervals themselves may be small. 

It is natural to think that a weight represents the average sensitivity of a coef- 
ficient of the nominal model with respect to the variation of data points. Thus, we 
establish the following reasonable algorithm for selecting weights. 

Suppose the test data consists of N data points obtained at corresponding fre- 
quencies, i.e., 

D{ju) := {D(ju/i) = a; +j(3i, i = 1,2, •••,#}. (3.3) 

Let us define the I th model set as follows: 


j D(jtui), i = l 

\ G^juJi), i = 1, 2, • • • , l - 1, l + 1, • • • , N 


(3.4) 


In other words, the model Gi(juj) is identical to the nominal identified model G T (juj) 
with the I th data point replaced by the Z th component of the test data D(juj). Now 
we construct the Z th identified model, which we call Gf(s ) which is identified from the 
I th data set Gi(jiv). Let 


g!(s) 


Too + + n l 2 s 2 + n l 3 s 3 + • • ■ + n l n s n 

d l 0 + d[s + d l 2 s 2 + d l 3 s 3 + • ■ ■ + d l n s n 


(3.5) 


and 



n Q Ui 



(36) 


If we assume that \Gi(juj) — Gi(ju>)\ is small, the sensitivity of the coefficients of the 
nominal model with respect to variations in the I th data point is described as 


dp 

dG T (juji) 


:= |n 0 -7io| 


\n r 


Til 


I do - 41 


\dn - 41 


(3.7) 


Collecting the sensitivity of the coefficients of the nominal model with respect to the 
variation of all the data points, Z = 1,2,*--, N, we have 


dp 

dG*(j ui) 


ap -i 
ao r (jw ,) 

9G'(ju, 2 ) 


§2 


K 41 

••• K-41 

••• 

14 

-41 1 


|n 0 - 4| 

K-41 


14 

-41 

(3.8) 

1 

o 

1 . 
3 

• • ■ K-41 

!*-<(? 1 

14 

-41. 
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The weights are then defined as the average of these for each coefficient: 


w := 


^[s£.l«o-»ol • 


E£,k-<l e£iK-<4I ••• 


W„o ••• W nn W^ w dn 


(3.9) 


3.3. Interval System Identification 


In this subsection, using the weight that represents sensitivity of the coefficients of the 
nominal transfer function, we develop an interval model that satisfies three conditions 
given above. After we determine an appropriate weight vector, we need to find £± 
to satisfy the given requirement. We now first consider the membership requirement. 
Recall the nominal system given in eq. (2.3) and substitute s = ju>, then we have 


G\ju) := 


n(juj) 
d(ju>) 
up + juni 
d 0 + jujdx 


u > 2 n 2 


co 2 d 2 


(n 0 — u> 2 n 2 + u> 4 ra 4 — •■•)+ j(wn x — + w a n 5 — 

(d 0 — uj 2 d 2 + u; 4 d 4 — ■■■)+ j(u)d x — u ; 3 d 3 + u 5 d s — 
n even (u>) + 7'n odd (u;) 


jui 3 n a + u; 4 ti 4 + juJ s n 5 
- ju 3 d 3 + u; 4 d 4 | j(j) 5 d a - 
) + j{unx — v 3 n 3 + w 5 n 5 

\ . •/ 1 3J , . .RJ 


n even (a>) + jn odd ( w) 
d even (u;) + jd odd (u;) ' 


(3.10) 


Since the nominal model transfer function G I (s) cannot perfectly represent the data 
set D(ju>), we have the following relationships for a particular frequency Wi. 


D(jWi) = OCi + jfii 

« G t (M) 

+ jn odd (LJi) 


d evcn (u>i) + jd odd (u>i ) ' 

The difference may be added to the coefficients of the nominal model as follows: 


(3.11) 


D(M) 


OLi + jfii 
(h 0 - u) 2 fi 2 + 


■) + j{win x - o;fw 3 + • • •) 


(d 0 - w?d 2 + •••) + j( w id i ~ w ?d 3 H ) 


(3.12) 


where 


hi := ni + vf ni £m 
di := dj+w^e^, 


Vi. 


If we rewrite this in terms of a linear matrix equation, we have 

A(u>i, cti, /3i)We l = B(wi,ai,f3i ) 


(3.13) 


(3.14) 
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where 


Oti -fiiWi 

Pi OL^i 


A <*;,&) := 

-C(iU)f PiU)f ••• — 1 0 u>j 0 

-fiiU)? — a, O'? ■ ■ ■ 0 (jJi 0 —Ujf 


-U)7 


W := 


0 


w dl 


w 


«o 


W„ 


£* := e: 




<0 


and B{u iy OL U Pi) : = 


-a; (d 0 - u;?d 2 + ufd 4 ) + Pi{ui t di - w?d 3 + w\d$ 

+(n 0 — u>fn 2 + c*^ra 4 ) 

~Pi(d Q — a pd 2 + vfd 4 ) - afaidx - ufd 3 + wfd 5 

+( 0 ;^! - a }fn 3 + u>fn 5 ) 


a;,” 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


Here we assume that without loss of generality A(a>i, has full row rank. It 

is well known in linear algebra literatures [10] that the smallest norm solution can 
be computed by 

^ = X(0 r [A(.)A(.f]-'S(0. (3.19) 

After finding for all i = 1, 2, •, V, we determine the interval as follows: 

e n* : = n Y n (0> 4J e n k '■= max{0, e^} (3.20) 

e d k : = m ? n {0) e dJ e d k : = max{0, e' dk } (3.21) 

for all k. Clearly, the procedure guarantees three requirement given earlier. 


4. MODEL VALIDATION 

Model validation can be accomplished by comparing the frequency domain charac- 
teristics of the obtained interval model and the test data. We use polar plots and 
Bode plots to accomplish this task. The polar and Bode plots of the test data is 
assumed to be given. Therefore we here only discuss the problem of obtaining polar 
and Bode plots of a proposed model interval system. The polar plot of an interval 
system is obtained by determining the image set of the set of transfer functions at a 
fixed frequency uj and taking the envelope of such sets obtained by sweeping over uj 
ranging from 0 — > oo. This is done via the following theorems. The proofs of these 
theorems may be found in [11]. Before we stating them some notation is necessary. 

Let G(s) be an interval transfer function family: 

G(„) := jcW = | N( S ) € jV(.), Z>(.) € V(.) J (4.1) 
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where Af{s) and V(s) are families of interval polynomials. The magnitude and phase 
of a transfer function G(s) at a frequency u are defined as Hg{w) and respec- 

tively. It is natural to define the maximum and minimum values of 0*0 and 4>g(v) 
at each frequency u>. 


and 


Let 


Hg( w ) := 

inf Hg{u) 
GeG 


Ag(w) := 

sup fiG{u) 
GeG 

(42) 

£ G H := 



^g(") : = 

sup <f>G{u)- 

(4.3) 

GeG 


» := iK(s) 

1 | i = 1,2,3, 4}. 

(4.4) 


where Kjf(s) is the i th Kharitonov polynomial associated with the interval polynomial 
Af(s) and similarly, we define Kd{s) to be the set of Kharitonov polynomials associ- 
ated with 2?(s). We now introduce a special set of segments that joint appropriate 
pairs of Kharitonov polynomials. These segments are known as extremal segments. 


Sjs(s) := 

{aJT'M + (1 - A )Kt(s) | A 6 [0. 1], (i, j) £ {(1, 2), (1,3), (2, 4), (3, 4)} (4 5) 

and Sd(s) is defined similarly. 

Introduce the following transfer function sets: 

Gk := I N(s) 6 /Cy(s), D(s) G /Cd(s)| (4.6) 

g b := I (JV(.),D(.)) € (JV(s) x 2 >(»))e} (4.7) 

where 

(JV(.) x DM). := {(JVM x DM) | 

N(s) 6 JC^M.DM € S„(s) or JVM € «#•(»), 7>M € £„(*)}• ( 4 '«) 


Theorem 1. [11] For every frequency u; > 0, 

Pg(v) = fiG B (v) 

£ G (“) = Hg b W 

Theorem 2. [11] For every frequency u> > 0, 

^g(w) = ^GkC^) 
= & Gk ( u ) 


8 



Theorem 3. [11] The image sei of the interval transfer function G($) evaluated at 
s = jw in the complex plane is bounded by the image of its corresponding extremal 
subsets: 

dG(ju) C G E (ju>) 

where <9(*) denotes the boundary of a set . 

Once the image set is generated at a fixed value of w — omegai one can verify if the 
given data point D[u>i) is contained in G (juii) or not. This verification can be carried 
out for every value of u;* to verify whether or not the proposed interval model explains 
the experimental data. 

In the next section, for the sake of illustration we apply these ideas to a large 
space structure experimental facility developed at NASA’s Langley Research Center. 

5. APPLICATION TO A MINI - MAST SYSTEM 

5.1. Model Description 

The Mini-Mast system shown in Figure 2 is a 20.16 meter long deployable truss 
located in the Structural Dynamics Research Laboratory at NASA Langley Research 
Center. It is used as a ground test article for the research in the areas of structural 
analysis, system identification, and control of large space structures. The Mini-Mast 
was constructed from graphite-epoxy tubes and titanium joints, by using precision 
fabrication techniques. The 102 measurements shown in Figure 2 were derived using 
51 noncontacting displacement sensors distributed from Bay 2 through Bay 18. Three 
shakers are located circumferentially around the truss at Bay 9 and their locations are 
selected primarily to excite the low frequency modes below 10 Hz. There are three 
modes which consist of two bending and one torsion mode in this low frequency range. 
These three modes are designed to be separated from the other frequency modes. 
The experimental data used in this example are obtained by using one displacement 
sensor output at Bay 9 from one input. In this example, we use the experimental data 
within the 45 radian/sec low frequency range with 180 frequency data points. This 
low frequency range covers the three low frequency modes described earlier. Figure 
3 shows the frequency response test data. 

5.2. Interval Model Identification 

Using the weighted least square method described earlier, we select W T (ju;) shown in 
Figure 4. The identified model obtained is 

qI,\ _ U 0 + U l 3 + + 7*5 S 5 

d 0 + dis + d 2 s 2 + d 3 s 3 + d 4 s 4 + d$s 5 + s 6 


where 
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no 

n i 
n 2 
n 3 

714 

n 5 

The eigenvalues of the identified model transfer function are as follows: 

Model no. Eigenvalues Mode 

1 —7.11 x 10 -2 ± y’5.356 first bending mode 

2 — 4.22 x 10 -1 ± j’26. 302 first torsion mode 

3 — 7.99 x 10 -1 ± jZ8 . 616 second bending mode 

The magnitude and phase comparison of the test data and the identified model 
are given in Figures 5 and 6. In Figures 5 and 6, the dashed lines denote the frequency 
response of D(ju>) and the solid lines denote the frequency response of G J (jw). The 
dotted lines in Figure 6 indicates the error in magnitude (i.e., \D(jw) — (j (yu>)|) for 
illustration. 

We now create intervals around this nominal identified model. The weight selec- 
tion method described in Section 3.2. gives the following weights for each coefficient: 

w no = 2.7053 x 10 w do = 3.9152 x 10 3 

w ni = 1.2041 w dl = 2.2715 x 10 2 

w n2 = 2.3214 x 10” 1 w d2 = 5.8095 x 10 

w n3 — 4.0113 x 10 -3 w dl = 1.5250 

w n „ = 2.4768 x 10- 4 w d4 = 5.9161 x 10~ 2 

w„ 5 = 2.9620 x 10“ 6 w d5 = 1.2520 x 10~ 3 

This set of weights produced the following interval system: 

_ rip + nis + n 2 s 2 + n 3 s 3 + ra 4 + n, 5 s 5 
d 0 + dis -f d 2 s 2 + d 3 s a + d 4 + d 5 s 5 + s 6 

where 

h 0 G [-6.0061,-5.6661] x 10 4 
h\ G [ 4.7696, 7.8401] x 10 2 
h 2 G [-8.9589,-8.6096] x 10 2 
h 3 G [-4.9791, 5.6904] x 10" 1 
n 4 G [—9.7491, —9.3729] x lO” 1 
h 5 G [-0.8570, 1.1648] x 10" 4 

5.3. Model Validation 

The following figures show that the interval model obtained here is a valid interval 
model for the given test data set. First, Figure 3 shows the polar plot of the test 
data for every measured frequency. Each mode of the polar plot has been separated in 


d 0 G [2.9361,2.9795] x 10 7 
d x G [2.0663,2.2335] x 10 5 
d 2 G [1.0923,1.0986] x 10 6 
d 3 G [2.6720, 2.8406] x 10 3 
d A G [2.2090, 2.2186] x 10 3 
d 5 € [2.5255,2.6190] x 10° 


= -5.78 x 10 4 d 0 = 2.96 x 10 7 

= 5.88 x 10 2 d x = 2.15 x 10 5 

= -8.74 x 10 2 d 2 = 1.10 x 10 6 

= 0.073 d 3 = 2.75 x 10 3 

= -0.967 d 4 = 2.21 x 10 3 

= 3.48 x 10" 5 d 5 = 2.58 
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Figures 7,8, and 9 for illustration. These figures show that every data point of the test 
data is bounded by the image set generated by the interval model at the corresponding 
frequency. Figure 10 was drawn for the entire frequency range. Similarly, Figures 11 
and 12 show the magnitude and phase plots of the test data and the interval model. 
Clearly, both magnitude and phase plots of the test data are contained in the tightly 
bounded tubes representing the boundary of the frequency responses of the interval 
system. 

6. CONCLUDING REMARKS 

A new algorithm to construct an interval transfer function from available frequency 
domain data of the plant is presented. This model captures the complete frequency 
domain characteristics that a single identified model cannot capture. This interval 
model can be directly useful to analyse the robustness properties of any proposed con- 
troller using the well developed theory of parametric robust control. Such an interval 
model can predict the worst case stability margins associated with the controller. 
This serves as a lower bound on the worst case stability margin of the actual (rather 
than the model) system. 
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